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ABSTRACT 



An analysis of the use of planing data to predict impact 
loads of a dihedral hydrofoil suitable for seaplane installation 
is presented and compared with experimental model tests. 

The experimental data verifies the theory very well for the 
prediction of maximum loads. No empirical correction factors are 
used and the theory agrees with experiment to an average error of 
6 , 2 % for the U7 runs analyzed. However, experimental time-to- 
peak tends to lag the theoretical, and the experimental draft 
time-history diverges from the computed after the point of maxi- 
mum load. 

One initial flight-path angle is considered for a comparison 
between the maximum loads predicted by these tests for a hydrofoil 
equipped seaplane and the maximum load expected of a chines-dry 
landing hull. 

One initial flight-path angle is considered to compare the 
relative efficiency of the two models tested as energy absorbing 
devices . 

In addition, a sample digital computer program useful in the 
application of the theory to the prediction of impact time- 
histories is presented. 

This investigation, conducted at the Experimental Towing 
Tank, Stevens Institute of Technology, Hoboken, N.J., is the 
thesis submitted by the authors in partial fulfillment of the 
requirements for the degree of Master of Science. 



NOTATION 



Mass of impacting body 

Virtual mass of displaced fluid 

Initial vertical velocity 

Instantaneous vertical velocity 

Vertical acceleration 

Mass density of water, 1.938 slugs/ft. 

Wetted half -width 
Length of float 

Wedge, dihedral or deadrise angle 

Virtual mass of flow plane at step per unit distance 
in keel direction 

Geometric angle of trim 

An area of fluid in a slice taken perpendicular to 
the keel 

Displacement in a direction normal to the keel 
Velocity normal to the keel 
Acceleration normal to the keel 

Displacement of the step from origin in a direction 
perpendicular to the keel 

Displacement in a direction parallel to the keel 
Velocity parallel to the keel 
Acceleration parallel to the keel 

Displacement of the step from the origin in a direction 
parallel to the keel 



Draft in y direction 

Draft of step in y direction 

Elemental force normal to the keel 

Force normal to the keel 

Total mass affected by passage of body 

/• 

M w = p cot T \ n w df 



Displacement in a vertical direction 



% 

M 

Force in vertical direction 
2 

COS T 

Horizontal velocity 
An arbitrary increment of time 
Initial flight-path angle 



Variation of the velocity potential with time in the 
general non-steady Bernoulli equation. 



A Study of the Hydrodynamic Impact Loads on 



a Dihedral Hydrofoil Suitable for Seaplane 

Use 



INTRODUCTION 



The hydrofoil as a seaplane alighting device has received 
serious consideration as a possible means of relieving rough- 
water and high speed landing operational problems. 

Compared to the conventional hull, the small area of the 
hydrofoil offers a considerable reduction of landing impact 
forces, while the fixed, relatively short hydrofoil chord limits 
center-of-pressure travel with a consequent reduction of undesir- 
able pitching moments. To a large extent these particular advant- 
ages of the hydrofoil are shared with the hydroski, but the hy- 
drofoil offers the additional advantage of lift/drag ratios 
superior to the ski, especially at low speeds . 

Although the hydrofoil offers much promise for seaplane use, 
with the exception of a single experimental report of a classi- 
fied proprietary nature, no studies of a general nature concern- 
ing its impacting properties are known. 

The object of this work is to provide a rational foundation 
for the stucfy and design application of impacting hydrofoils for 
seaplanes. To accomplish this end, an examination will be made 
of the use of planing data to compute impact loads and the 



practicability of extending its use to the specific case of the 
hydrofoil; then to analyze and compare such analytical results 
with experimental model tests. 

This investigation has been conducted by Lieutenents Richard 
A. Hoffman and Arthur C. Derrick, U.S. Navy, Naval Postgraduate 
School students at the Experimental Towing Tank, Stevens Institute 
of Technology, Hoboken, New Jersey. Special acknowledgement is 
due to Professor B.V. Korvin-Kroukovsky, Mr. Daniel Savitsky, 

Mr. P. Ward Brown and Mr. Robert Van Dyck of the Experimental 
Towing Tank Staff for their most helpful interest and guidance in 
this work. Thanks are also in order to all of the other members 
of the staff who so freely gave of their assistance, experience 
and advice during this investigation. 



EXAMINATION OF THE USE OF PLANING 



DATA TO COMPUTE HYDROFOIL IMPACT LOADS 



Historical Background 



The earliest impact theories for the landing of seaplanes 
were based on the work of von Karman (Ref. l) in which he con- 
sidered a wedge shaped body dropped vertically into a horizontal 
water surface. With the assumption that momentum was transferred 
to a hypothetical associated mass of water attached to the sea- 
plane and that the total momentum of the system remained constant, 
von Karman was able to write the momentum equation: 



Mz = M z + m z 
o w 

where M - mass of the impacting body 

m w - virtual mass of displaced fluid 

o 

z q - initial vertical velocity 
z - instantaneous vertical velocity 
With the further assumption of two dimensional flow in planes 



normal to the keel, von Karman considered at each instant the 
virtual mass is equal to the mass of water contained in a semi- 
cylinder of length equal to the length of the float and of 
diameter equal to the width of the wedge at the water surface 
(Fig. l). This expression for the virtual mass substituted 
into equation (l) led to the force equation: 



(1 



a. 



F 




( 2 ) 



where c - wetted half width 
JL - length of body 

Subsequently Wagner (Ref. 2) took into account the wave 
rise at the wedge surface and defined the diameter of the semi- 
cylinder of virtual mass on the basis of deadrise and draft 
(Fig. la)j and Pabst (Ref. 3) introduced an aspect ratio factor 
to account for end losses in approximating oblique impact on the 
basis of the two-dimensional calculation. 

Mayo (Ref. U) in 19h% analyzed the basic von Karman theory 
and concluded it erred in assuming the total momentum of the 
system remained constant in that it did not account for momentum 
transfer onto the wake. Mayo pointed out that a more realistic 
force equation would be 



in which the rate at which momentum is imparted to the downwash 



F 



y 



d / • s • • 

-rr (m y) + m x y 
dt w J s J 



(3) 



where y - velocity normal to keel 

x - velocity parallel to keel 



m - virtual mass of flow plane at step per unit dis- 
tance in keel direction 



is included 



Mayo, using equation (3), was able to apply the von Karman 
theory to chines dry oblique impact by imagining the fluid be- 
neath the wedge divided into slices of length dx by planes 
drawn normal to the keel. As the wedge moves through these 
slices, the flow in each slice is assumed to be identical with 
the flow in the case of the wedge dropped vertically. By assum- 
ing the fluid in each slice to be independent of the motion in 
adjacent slices, the total force on the wedge is equal to the 
sum of the forces of all the slices acting on the wedge plus the 
rate at which momentum is imparted to the downwash by flow planes 
sliding off at the step. 

Mayo pointed out that in planing the first term of equation 
(3) is 0 and the hydrodynamic force is equal to the last term 
and therefore planing data might be used to evaluate the momentum 
of the flow in the normal cross sections. 

Subsequently, in 19U8 Steiner (Ref. 3) followed Mayo's ap- 
proach to obtain correlation between planing data and one of the 
virtual mass expressions, using limited data mostly confined to 
the chines dry case. In 19!?1 Smiley (Ref. 6) used planing to 
obtain a correlation between planing and impact for a flat plate, 
obtaining an impact pressure coefficient in terms of an "equiva- 
lent planing velocity". 

With the advent of heavily loaded high length-beam ratio 
hulls which immersed chines in landing and of flat bottom skis, 
certain inherent weaknesses of the classic theory precluded its 



use for these modern applications. For example: for a flat 

bottom the force equation (2) yields an infinite force; while 
for the case of chines wet impact, the virtual mass becomes 
constant at chine immersion and does not increase with increasing 
draft. 

In 19!>2 Schnitzer (Ref. 7) treated the case of chines wet 
impact by constructing a composite deflected mass from the de- 
flected (virtual) mass apparent at chine immersion plus a 
deflected mass derived from Bobyleff's theory of the flow about 
a two-dimensional V-shape of finite width. However, Schnitzer* s 
theory was lacking somewhat in generality in that it attempted to 
absolutely define the geometry of the virtual mass. 

Also in 1952 Smiley (Ref. 8) presented an analysis for the 
landing impact in which the landing loads and motions are ex- 
pressed in terms of the planing properties of the body. However, 
Smiley reverted to an empirico-theoretical virtual mass expression 
at chine immersion. 

In 195U, Brown (Ref. 9) exploited to the full the concept 
that impact characteristics may be predicted with the use of 
planing data, avoiding the unverifiable definition of the geome- 
try of the virtual mass. 

Brown's approach is especially attractive for the computa- 
tion of impact loads on such complicated shapes as hydrofoils and 
skis, since it is unhampered by empirico-theoretical expressions 
for virtual mass. In its published form however, the derivation 



of Brown's theory breaks down when the leading edge of the body 
becomes wetted as in an impacting hydrofoil. During the course 
of this work Mr. Brown became associated with the Experimental 
Towing Tank and through his cooperation the derivation was modi- 
fied so that the theory is applicable to any impacting body for 
which planing data is known. 

Because of its logical simplicity and its applicability to 
any planing body. Brown's approach is the one followed in this 
analysis and the revised derivation is discussed below. 

Seaplane Impact Theory as Applied to the Hydrofoil 

The following analysis considers the impact of a straight- 
keeled planing body of arbitrary cross section. 

The following assumptions are made: 

1. Trim, t , remains constant throughout the impact. 

2. Inertia forces predominate. 

3. The body is a weightless mass, weight supported by wing 



U. The virtual mass M w , is some function of the draft. 
Consider the impacting body shown in Fig. 2. In the fixed 



lift. 



plane AA' a cross sectional area of water n is influenced by 

w J 

the passage of the body. 



The momentum of the fluid in the 



element n normal to the keel 
w 




(U) 



The area n is considered to be a function of the draft f , 
w 

since in a booty of known geometry at fixed trim, ary cross- 
sectional area may be determined if the draft is known. When the 
body has passed through AA' , the momentum in the element is con- 
sidered as shed to the wake. This concept is discussed more fully 
in Monaghan (Ref. 10). 

The rate of change of momentum ^ 
in an element fixed in space J 



Since 



at (p n „ Jx) 



dn 

- eJ*<vs, ♦ y s -dT> (5) 



dn 



but 



w 



. djr 

djf dt 



dt 

and from Fig. 2 J = y - x tan t 

s 

and since the plane is fixed x is constant 

* f = 7 

. . j s j g 

2 dn 

and (£) becomes = p £x (y g n w + y g -gy) 

Rate of change of C Force exerted on 

momentum in an y = a body by element 
element 



(6) 



o dn 

V ■ p 5x <?. \ + d#> 



( 7 ) 



This determines the elemental force acting on the body in a 



fixed plane. To determine the total force acting on the body at 
any one instant of time, consider the body's position, velocity 
and acceleration to be instantaneously fixed. Now let the 



9 



previously fixed plane move along the body, recalling as x 

varies, S varies and with J 3 , n^ „ 

Since J 3 = y - x tan t , under the above conditions 
s 

d = - tan t dx and 



x +L 

„ f S / , * 2 w\ , 

F y * ( p n w ♦ p a- j ] ** 



f 

2 r s dn 

= p y s cot t \ n w d f + p y s cot x ^ jy d f 



( 8 ) 



Now define M^ , the total mass effected by the passage of 



the body as: 



f s 

M = p cot t ( n df 
w K \ w J 



dM fs dn 

= P cot x ■— df 

J o 

Equation (9) may not be readily apparent, but it may be proved 
by a series expansion of n as a function of JT . 

Therefore, substituting (9) into (8) 



(9) 




M 

w 




dM 

w 

d 4 



From Figo 2, f ^ = z sec t 

df = sec t dz 



( 10 ) 



and equation (10) becomes 




M 



w 




dM 

w 

dz 



COS T 



( 11 ) 



From Newton’s Second Law 




»v 



M 



w 



dM 

w 

dz 



cos T 



y= - 



o dM 
• 2 w 

V -r— COS T 

•’s dz 



M 

M(1 + f ) 



M dM 

Define , then — ~ = M 

M *• ’ dz dz 

Substituting in (12) 



•• ° 2 d-K. _ 1 

y s = y s to cos T TTT57 



( 12 ) 



(13) 



df i 

Determination of dz and 

Consider planing as a special case of impact, Fig. 3» In 
this case 



F = F cos t 

2 y 



y g = V H sin t 



where is the horizontal planing 

velocity 



and substituting in (ll) yields 



11 



F = V 2 sin 2 x cos 2 x M 
z H dz 



= z_S_ 



dz V 2 sin 2 x cos 2 x A 
H o 



= 



_g_ 



2 2 2 1 
V„ sin a cos x A _/ 
n O O 



F dz 
z 



(Hi) 

(15) 



The constant of integration is zero since ^ = 0 at z = 0 

Thus the virtual mass ratio is defined in terms of the 

integrated lift/draft curve and constants. 

In order to change the equation of motion (13) which is in 
terms of velocities and accelerations normal to the keel to the 
more useful terms of vertical velocities and accelerations (as 
are measured in the towing tank), substitute in (ll) 



F cos x = F 
7 z 

y = V„ sin x + z cos x 

S n 



y = z cos x 



dM 

• •p . • \ 2 -rr 0 

* F = z cos x M + (Vu sin x + z cos x) — — cos x 
• * ?. "w rl ^ ** 



dz 



( 16 ) 



Since F = - M- , combining with ( 16 ) and (lii) and 
z z 

rearranging yields 



- z = F_ 



Vjj sin x + z cos x 



Vy sin x 



M(1 + ^ cos t) 



(17) 



t 



2 

= cos t yields 



Letting ^ 



z = F 



sin t + z cos t 
V H sin t 



_1 

M (1 



(18) 



The resulting non-linear equation ( 18 ) does not lend itself 
to a closed solution, therefore a step-by-step integration is 
necessary to compute the time history of impact. 



Step-by-Step Integration of Equation of Motion (18) 

The purpose of this integration is to determine the time- 
history of hydrofoil impact after striking the water under known 
initial conditions. During this impact the trim is assumed to 
remain constant, as is the horizontal velocity. 

The horizontal velocity, trim angle, initial vertical veloc- 
ity, weight of the body and a curve of steady state planing lift 
versus draft for the horizontal velocity of impact are known. 

Then to perform the integration the only quantity needed is a 
time interval (At) small in respect to the time from contact to 
maximum load. 

If zero time is the instant of contact, the average draft 
during the first interval is 

z At 

d = ° 0 - - 

av^ 2 

Then from a curve of planing lift versus draft for the foil, 

the lift (F in equation 10) for the average draft is determined 
z 



and from an integrated plot of F versus draft, the quantity 

z 
-z 



£ 



F dz used in the determination of u' is found, 
z 



To compute the acceleration, the planing lift is multiplied 
by the quantity 



& 

A. 



V„ sin t + z cos t 
H o 



x sm t 



1 + 



g 



2 2 

A V„ sin 4 t 
o H 



£ 



F dz 
z 



all quantities of which are now known. 

The vertical velocity at the end of the interval is 



z., = z - z At 
1 o 



and the draft is 



z + z, 

i o 1 „ , 

d l a - “ 



Using z^ and d^ as initial conditions, the average draft 

for the next interval would be 

z.. At 

v, ■ d i * - v 

and the procedure is repeated. 

The most efficient manner to accomplish the integration is 
in tabular form, a sample of which is given as Table I. However, 
even at best the step-by-step integration is a tedious and time- 
consuming task, and since each result is dependent upon the 
previous one, it is a process in which errors accumulate. 



An examination of Table I reveals that the step-by-step 
integration is in a form to be easily programmed for a digital 
computer. This was done for the ELECOM 100, and a sample flow 
chart easily adaptable to any digital machine is included as 
Table II. 

Computation of Impact by Digital Computer 

The sample flow chart is so constructed as to be useful for 

the determination of impact time histories of any arbitrary bodies 

for which planing data is known. 

The main program is standard and need not be altered if the 

lifting characteristics of the body are changed. The F vs. z 

z 

r z 

and \ F z dz vs. z functions are in subroutines which may be 

changed with comparative ease. 

With this sample program the only inputs needed in addition 

to the above subroutines are V , z , sin t , cos t , A , g , 

no O 

At , 1 and 2. Therefore numerous runs may be computed in a very 
short time . 

All theoretical curves except for the run of Table I were 
computed using this program. 

EXPERIMENTAL INVESTIGATION 

Purpose 

While the planing analogy as a means of computing impact 
loads has become well established in recent years, Refs. £ through 



9 and 11,, all known experimental data analyzed by this technique 
has concerned hulls or skis of low aspect ratio and a dry leading 
edge. Although intuitively it might be expected that the analogy 
would hold for any planing body for which the lift coefficients 
are known, the large aspect ratio and submerged condition of the 
hydrofoil make it so different from the bodies previously studied 
as to indicate a definite need for experimental data to verify 
the applicability of the use of planing data to compute impact 
loads. To this end the facilities of the seaplane tank. Experi- 
mental Towirg Tank, Stevens Institute of Technology were employed. 

Equipment and Procedure 

Two hydrofoil models made from half-hard brass bar stock 
were used for the tests. The first model consists of a center- 
supported V-shaped hydrofoil with a constant chord of one inch 
and a length of each arm of five inches. The dihedral of the 
arms is 30 degrees from the horizontal. The bottom of the foil 
is flat with sharp leading and trailing edges. The center sup- 
port is l£ inches long and sharp-edged to prevent excessive drag 
from developing as a result of spray. Bracing from the foil tips 
to the center strut were added to insure rigidity. See Fig. U. 

The second hydrofoil is the same except that the chord, in- 
stead of being constant, has a taper ratio of 1:2 with a one inch 
chord at the center and a two inch chord at the ends. See Fig. U. 

The tests utilized the Tank 3 "Nutcracker'’ apparatus with a 



Schaevitz type AB-2 accelerometer with "damped natural frequency" 
on the order of 90 cycles per second and approximately 0«7 criti- 
cal damping. The standard Tank 3 speed and recording equipment 
was used in connection with two Hathaway oscillograph recorders 
for time -variable measurements. 

The readings of heave and acceleration are picked up by and 
transmitted by means of Schaevitz electro-mechanical differential 
induction transformer units which give approximately linear re- 
sponse of output voltages to input mechanical deflections. 

The Tank 3 "Nutcracker" apparatus (Figs. £, 6) consists of a 
pole free to move vertically through two collars housing low- 
friction ball bearings but limited from horizontal rotation by 
means of a pivoted arm and elbow linkage. This apparatus has a 
built-in Schaevitz unit to measure heave of the pole by means of 
a string wound around a pulley on a worm gear. A solenoid- 
controlled heave lock device is built onto the "Nutcracker". In 
the tests this device was tripped electrically as soon as the 
carriage reached a predetermined steady speed energizing the 
solenoid which pulled a pin out of the pole and allowed the model 
to drop into the water. The "Nutcracker" was originally fitted 
with a spring and cam unloader which provided a constant lift 
force independent of the model's vertical position without add- 
ing any appreciable mass to the model. 

However, after a number of trial runs using the spring and 
can unloader, examination of the oscillograph tape revealed that 



the unloader induced such an extraneous vibration into the system 
as to make analysis of the acceleration time-history impossible. 
Therefore, the weight of the model, pole and appendages was un- 
loaded by "negator" type sheave and line constant force springs 
during the entire range of heave, and downward velocity was 
imparted by using a ten pound weight on a pan attached to the 
pole by a "C" clamp. The weight was attached by nylon cord to 
the fixed (in heave) carriage, Fig. 7, so arranged as to be 
lifted from the pan at a fixed point before the model hit the 
water. See Fig. 8. Thusly the model was given downward momentum 
by the weight and when the weight was lifted off the pan by the 
cording, the model, being unloaded to a zero "g " condition by the 
"negator" springs, continued downward at a constant velocity 
until it struck the water. Examination of subsequent oscillo- 
graph tapes revealed that this system, although hurridly impro- 
vised due to the limited tank time available, did indeed impart 
a constant vertical velocity after the weight was lifted off and 
no discernible vibration to the accelerometer. 

The seaplane tank. Tank 3, where the tests were conducted is 
313 feet long, 12 feet wide and 6 feet deep when the water level 
is just even with the opening in the standpipe inside the tank. 

A maximum carriage speed of 50 feet per second can be attained 
with the selected speed, amplidyne controlled, carriage driving 
motor. A tachometer generator, the armature of which is driven 
by the carriage driving motor shaft, supplies electrical information 



to the field of the amplidyne system to provide speed control. 

With a selected voltage, representing a certain speed, the 
carriage speed can be maintained to within 0.06 feet per second. 

The speed of the towing carriage is measured and recorded by 
an electronic tachometer specially designed and built by the 
General Electric Company. The accuracy of this instrument is in 
the order of a quarter of one percent of the speed being measured, 
in line with the degree of accuracy required in the drive system. 

The instrument receives its impulses for speed measurement 
from photoelectric scanning discs mounted on the cable-drum shaft. 
By suitable selection of the number of scanning holes and use of 
band-spread switches, a small increment of speed can be selected 
and expanded to cover the full scale of the instrument so as to 
insure a high degree of accuracy of speed measurement. 

As an auxiliary means of determining carriage speed, a mech- 
anical "chronograph" device is used in which the distance travelled 
by the carriage during measured time intervals is plotted on mov- 
ing paper tape. 

The models were installed in the overhead towing carriage as 
shown in Figs. 5, 6, 7 and 8 with provisions made to set the angle 
of trim by a spirit level protractor. With the trim angle set at 
zero degrees, adjustments were made to insure that there was no 
angle of roll or yaw. By lowering the model into the water and 
marking the free surface intersection on either side, roll was 
considered zero when the wetted length of foil span on either 



side of the centerline was equal to the other. Yaw was checked 
by attaching a straight machined bar with its lengthwise axis in 
the direction of motion to the model strut. Then, if by moving 
the carriage forward the edge of the bar remained at a constant 
distance from the tank wall, yaw was considered zero. 

Calibration in heave was made by moving the pole up and down 
in one inch increments. Heave tape readings were made at each 
position. Therefore inches of tape per inch of actual heave 
could be determined. 

Calibration in acceleration was made by detaching the accel- 
erometer from the pole and standing it upright (zero "g"), on its 
side (minus one "g") and upside-down (minus two "g"), and taking 
tape readings at each position. Thus inches tape per "g" could 
be determined. 

Calibration in heave and acceleration was made daily before 
starting runs and whenever a change in set-up was made. 

In addition, before each run heave was marked on each tape 
at the bottom stop; where the trailing edge tips of the model 
just touched the water; and where the trailing edge apex of the 
model just touched the water, from the intersection of the apex 
touching trace with the running heave trace, the beginning and 
end of impact was determined. 

Runs were made at a constant forward speed of about 
30.7 ft. /sec. and for each trim angle, various vertical veloci- 
ties, see Table III. The vertical velocities in each case could 



not be precomputed due to friction in the pole and the stretching 
of the nylon cord which varied the length of drop over which the 
weight acted. All vertical velocities were computed from the 
heave vs. time oscillograph tape. 

During the planning phase of these tests it was concluded 
that in full scale landings full ventilation of the upper foil 
surface would undoubtedly occur. In order to simulate actual ex- 
pected conditions, the experiments were planned to achieve separa- 
tion at the leading edge and consequent ventilation by a combina- 
tion of high forward speed and a sharp leading edge. 

As demonstrated by Wu (Ref. 12) with full ventilation the 

lift coefficient of a submerged planing flat plate may be pre- 

2 'TT* sin 'T 

dieted by the Rayleigh formula ^ + f7's±n % ’ which can 

be modified for the effect of dihedral on angle of attack to 

yield C T = i— ? - 003 After a few trial impact runs it 
J L U +7rsin t cos p * 

became evident that the modified Rayleigh lift coefficient was 
not being realized at low trim angles. It was concluded that 
this was due to the fact that the models, being designed for 
rigidity and economy of fabrication, had too large a leading edge 
angle to achieve separation at low trims, and that at trims below 
the critical separation trim the upper surface was actually gen- 
erating negative lift. Since funds were not available to recon- 
struct the models, the trim angle was increased to a minimum of 
12° for the tapered-chord foil (Foil #l) and 10° for the constant- 
chord foil (Foil #2), where positive lift was attained. It was 



further concluded separation and subsequent ventilation occurred 
at and above these trims and they were the minimum examined in 
the impact condition. 

In order to obtain accurate lift data for the impact condi- 
tions, planing data was obtained experimentally for both foils 
at the trim angles and horizontal speed considered in impact. 

This experimental lift vs. draft data is plotted in Figs. 9 and 
10 for Foils 1 and 2, respectively. This planing data was ob- 
tained by the same apparatus. After unloading the model and pole 
with "negator" springs to a zero "g" condition, known weights 
were attached to the model and it was run down the tank at the 
forward speed to be examined in impact in a free-to-heave condi- 
tion. When the predetermined steady speed was reached, the heave 
(model draft) was read from the visual scale of the oscillograph. 
The plots of lift vs. draft were obtained by varying the load. 

COMPARISON OF THEORETICAL RESULTS 
WITH EXPERIMENTAL DATA 

General Discussion In order to determine the applicability and 
degree of accuracy of the developed impact theory to the specific 
case of the dihedral hydrofoil, a comparison was made with ex- 
perimental model test data. The results of the tests are pre- 
sented as: (l) experimental curves of vertical steady-state 

planing force vs. draft, Figs. 9 and 10 for foils 1 and 2 respec- 
tively, made at the single horizontal velocity and at the various 



trim angles examined in impact. The trim angles used were high 
enough to insure ventilation in the planing condition; (2) curves 

r z 

of the function IF dz vs. draft obtained from graphical inte- 
J o 

gration of the lift/draft curves. These curves are presented as 
Figs. 11 and 12 for foils 1 and 2 respectively; (3) plots of maxi- 
mum loading per run vs. initial flight-path angle. Figs. 13 and 
111 for foils 1 and 2 respectively; (U) tine-history comparisons 
with results obtained from the stepwise integration of the equa- 
tion of motion, equation (18). See Figs. l£ through 32; and (£) 
loop curves to compare energy absorbing efficiency of the two 
models at one initial condition, see Fig. 33* 

Experimental Data The experimental data obtained in the tests 
were: (l) curves of steady-state planing force vs. draft necess- 

ary in the computation of the theoretical time-histories and 
(2) experimental time-histories of draft and vertical accelera- 
tion. Experimental draft and acceleration were measured directly 
from the oscillograph tape. The initial vertical velocity, meas- 
urement of which is extremely important for an accurate comparison 
between theory and experiment, was computed from the slope of the 
draft trace just before water contact. Correct interpretation of 
the slope of the draft trace is quite dependent upon the judgment 
of the observer. For this reason, the initial vertical velocities 
were computed independently by both investigators and the recon- 
ciled values were used in this report. 



Comparison with Theory Time-histories of experimental draft and 



acceleration for 18 of the runs listed in Table III are compared 
in Figs. 15> through 32 with those computed from equation (l8) for 
the same initial conditions. The runs compared were chosen at 
random and are considered representative of the three trims used 
and approximately cover the range of initial flight-path angle, 

, used in the tests. Again, it is emphasized that full 
ventilation is expected in full scale landings and that the test 
impact data is for the case of upper surface ventilation. 

Results and Discussion 

Acceleration Examination of Figs. 15 through 32 reveals that 
the presented theory is in excellent agreement with the experi- 
mental values of maximum acceleration. However, a general tend- 
ency of the experimental time -to -peak to lag the computed is 
apparent. Although this lag could be blamed on the difficulty- 
in establishing an accurate reference line for acceleration 
measurement, it is felt that this time lag is a manifestitation 
of the inherent weakness of the planing analogy, i.e., the treat- 
ment of impact as a quasi-static phenomenon. It is well known 
that pressures on submerged objects as determined by the general- 
ized Bernoulli equation: 

T“it +gZ + p" F(t) 

bj. 

b t 



contain a 



term which increases in importance with increasing 



rapidity of time change. For example, in the two-dimensional 



wedge impact treated by "expanding plate theory" of von Karman 

and Wagner, the major part of the force results from the ^ 

a t 

term. However, this theory, as most potential impact theories, 
treats the particular case of uniform penetration velocity in 
which case the term furnishes the major resisting force. 

a u 

In decelerated penetration, as in a seaplane landing, the 
term would be reduced by the inertia of the decelerating virtual 
water mass. In the case considered in the present report, the 
consideration of impact as quasi-static can be considered as hav- 
ing provided a satisfactory first order answer. With increasing 

impact velocity however, the effect of the 7 -^ term as affected 

o t 

by deceleration can be expected to become more and more important, 
manifesting itself in a lag of the resisting force. As penetra- 
tion velocity is reduced and approaches zero the conditions 
approach static planing and satisfactory agreement of maximum 
acceleration is expected and in fact obtained. 

It is felt a more precise and rigorous impact theory for 
ventilated hydrofoils can be obtained from potential flow theory. 
Green's paper. Ref. 13, provides potential equations for the 
discontinuous flow about a flat plate as does the work of 
Rayleigh, Ref. lli, and Milne-Thomson, Ref. l£. However, all of 
these references express the potential function in terms of an 
intermediate mapping plane and the mathematical difficulties in 
expressing the potential in terms of the physical plane precluded 



further examination of this approach in the time available. 

It is felt that further examination of the application of the 
time dependent potential flow theory to hydrofoil impact will be 
most fruitful. 

Draft Comparison of the computed and experimental draft curves 
of Figs. 15> through 32, show good agreement during initial pene- 
tration. However, all experimental values of maximum draft are 
larger than those computed. It is felt that part of the reason 
for this deeper penetration is due to the time-lag effects dis- 
cussed above. Due to the time-lag, the energy is absorbed more 
slowly than predicted, resulting in a deeper penetration. 

Another reason for the variation is considered to be due to 
separation and ventilation effects. As the vertical velocity 
becomes small near maximum draft, the flow tends to close over 
the upper surface of the foil, and for these particular foils 
this reduces the lift. With less lift to resist penetration, 
the draft is increased. 

The ventilation effects are especially evident at low trim 
angles. For foil 1 at 12° trim and for foil 2 at 10° trim, 

Figs. 1 $ through 17 and Figs. 2U through 26 respectively show 
experimental draft values diverging from the computed after 
maximum draft, while at higher trims the experimental curves 
generally parallel the computed after maximum draft. Since loss 
of proper ventilation is more critical at the lower trims, the 



divergence of the experimental draft from the computed low trim 
indicates definite closing in of the flow at low trim and probable 
closing in effects at the higher trims. 

It is felt that these pronounced ventilation effects are due 
in a large part to the poor hydrodynamic design of the model foils 
tested. The authors of this work have proposed tests of a dihedral 
hydrofoil constructed for a dynamic model of an actual aircraft 
to the Experimental Towing Tank, and the results obtained therefrom 
should be in much better agreement with theory. 



Variation of Maximum Loading with Initial Flight-Path Angle 
Figs. 13 and lh reveal that the maximum loading of a dihedral 
hydrofoil increases linearly with initial flight-path angle, and 
with a slight increase of maximum loading with increasing trim. 



Effect of the Virtual Mass Term 



( i + V ) 



of Equation (18) 



Analysis of ten hand calculated step-by-step integrations of the 

equation of motion, of which Table I is a sample, reveals that 

the — 7 — r term of equation (18) is conservative, with a mini- 
1 +-*{_ 

mum observed value at maximum loading of .90 . Therefore its 
omission in initial estimate of maximum loading, which greatly 
simplifies hand calculation, should result in a maximum over- 
estimate of maximum loading of approximately 10$. 



Usefulness of Sample Computer Program of Table II Comparison of 
results obtained from the digital computer program of Table II 
with hand calculated results revealed the use of the computer 



resulted in a much greater degree of accuracy than possible with 
a slide rule at a great saving in time. Calculation of the com- 
plete time-history of impact by hand took approximately two hours 
Calculation by the digital computer took about 3 minutes of com- 
puter running time. 

Possible Application Although the hydrofoil models of these 
tests were designed primarily for verification of the use of 
planing data to compute hydrofoil impact loads, the general con- 
figuration of the models is envisioned as being suitable for 
seaplane installation. To afford some first-order degree of 
approximation as to the advantages of such a hydrofoil system, 
consider the models to be scaled up by a factor of 11.9. This 
corresponds to an aircraft having a landing speed of 63 knots 
and a gross weight of 9050 lbs. for foil 1 and 8000 lbs. for 
foil 2. Since these weights and speed are comparable to the 
Grumman JRF, the $ foot beam and 30° deadrise angle of the JRF 
were selected as further comparison parameters. Using Milwitsky' 
tables for chines-dry impact, Ref. 16, and the 10.3° flight-path 
angle of runs 1 and 37 for foils 1 and 2 respectively, maximum 
load factors for the chines-dry hull impact were computed as 
7.3"g" and 7.6 M g" respectively. Comparison of these predicted 
loads with the hydrofoil loads as listed in Table III for the 
same initial conditions show that foil 1 provides a X£>% and 
foil 2 a %6% reduction in maximum impact load over the chines-dry 



hull. This reduction is accomplished by a full scale foil of 
8.5 ft. span and the following chord dimensions: 

For foil 1 - an apex chord of 11.9 inches and a tip chord of 
23.8 inches. 

For foil 2 - a constant chord of 11.9 inches. 

To consider the foils as energy absorbing devices, the loop 
charts of Fig. 33 were constructed. Graphically integrating the 
area under the curves to the point of maximum penetration checked 
the initial kinetic energy, M z q ^ , within $%, Furthermore, 
comparing the area under the curves with the square area of an 
ideal shock absorber showed both foils to have an ideal effici- 
ency of approximately $6%, However, taking the ratio of rebound 
energy (area under lower limb of loop curve) to input energy (area 
under upper limb), showed that foil 1 returned only 17.5^of the 
original kinetic energy to the system as rebound while foil 2 
returned 28 . h%» Therefore, a tapered foil corresponding to the 
dimensions of foil 1 would provide a more damped landing device 
with lesser likelihood of the heavy damage experienced with the 
rebound impacts of the conventional hull. 

Conclusions An analysis of the experimental data obtained dur- 
ing hydrodynamic impact of a dihedral hydrofoil resulted in the 
following conclusions: 

(l) In general, the use of the quasi-static planing analogy 
to compute hydrofoil impact loads gives a good all-around approach. 



In particular, computed maximum loadings are in very good agree- 
ment with experiment, showing an average error of 6 , 2 % for 47 
runs. 

(2) Experimental time -to -peak tends to lag the computed. 

(3) Correspondence between computed and experimental draft 
is good until maximum force, after which a divergence is noted 
due to the time-lag of (2) above and separation effects. 

(4) The maximum impact loading of a dihedral hydrofoil 
increases linearly with initial flight-path angle and increases 
slightly with increasing trim. 

(f>) The effect of the Virtual Mass Term ( ^ J of equa- 
tion (18) is conservative, with a minimum observed value at 
maximum loading of .90. 

(6) The sample digital computer program of Table II improves 
the accuracy of the calculation of draft and acceleration time- 
histories and reduces the time required for such calculation 
from 2 hours to 3 minutes. 

(7) Compared with chines-dry impact at a 10.3° flight-path 
angle, the dihedral hydrofoil system offers a reduction of maxi- 
mum impact loading in the order of 5 0 %, 

(8) At approximately 8.9° initial flight-path angle, both 
foils tested had an ideal shock absorber efficiency of $ 6 % ; how- 
ever the larger tapered foil 1 returned but 17 , 6 % of the original 
kinetic energy to the system as rebound energy, while foil 2 re- 
turned 28.42, and therefore a foil similar to foil 1 is recom- 
mended for its damping qualities. 
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Sample Flow 'Chart for Step-by-Step Computation of Equation 
of Motion by use of a Digital Computer 

* Subroutine Shown is for Foil 2 at t = 12° 
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TABLE II 
Part B 

Cell Fill and Working Storage for Sample Digital Computer Program 
to Determine Time History of Impact 

* Main Program fill Subroutine Fill 
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TABLE III 

EXPERIMENTAL IMPACT RUNS 
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